Atmospheric, long baseline, and reactor neutrino data constraints on 6*13 
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An atmospheric neutrino oscillation tool that uses full three-neutrino oscillation probabilities and 
a full three-neutrino treatment of the MSW effect, together with an analysis of the K2K, MINOS, 
and CHOOZ data, is used to examine the bounds on #13 . The recent, more finely binned, Super-K 
atmospheric data is employed. For L/E^ > 10 4 km/GeV, we previously found significant linear in 
S13 terms. This analysis finds #13 bounded from above by the atmospheric data while bounded from 
below by CHOOZ. The origin of this result arises from data in the previously mentioned very long 
baseline region; here, matter effects conspire with terms linear in #13 to produce asymmetric bounds 
on 0i3. Assuming CP conservation, we find #13 = — 0.07±g;if (90% C.L.). 
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The phenomenon of neutrino oscillations [3, 0, 0, 0] has 
been observed in a variety of experiments: solar, long 
baseline (LBL) reactor, atmospheric, and LBL acceler- 
ator experiments. Including the constraint imposed by 
the CHOOZ reactor experiment one may quantita- 
tively determine the three mixing angles and two mass- 
squared differences that parameterize three-neutrino phe- 
nomenology Q. An outstanding question is the value 
of the mixing angle #13. Present analyses @| yield 
1 013 1 5= 0.15. We examine the impact of small effects, 
particularly those linear in 13 0, H, 0, [13] > on extract- 
ing this small parameter from the data. We find that 
including the full three-neutrino oscillation probabilities 
and a full three-neutrino MSW calculation are important 
for determining this mixing angle. 

Knowledge of #13 is a particularly important part of 
neutrino oscillation phenomenology because its value sets 
the magnitude of possible CP violating effects and the 
size of effects that might be used to determine the neu- 
trino mass hierarchy. There are presently three new re- 
actor experiments under development which are designed 
to measure #13, Daya Bay 11 1, Double CHOOZ [13], and 
RENOjil, as well as two long baseline experiments, 
T2K [lj and NOvA [15j . The subsequent generation of 
experiments, e.g., those which will ascertain the level of 
CP violation, cannot proceed until the current generation 
better determines the value of 6*13. 

The standard model of neutrinos conserves flavor, as 
is required by the data. Neutrino oscillations require 
adding a posteriori a mass term to the standard model 
Lagrangian. The standard model Lagrangian is diago- 
nal in flavor; the added mass term is diagonal in the 
basis which governs vacuum propagation. The relation 
between the two bases is given by a phenomenological 
unitary matrix U a i. In the absence of CP violation, it is 
real. We employ the standard representation [l(| written 
in terms of the three mixing angles, 6*12, 6*23, and #13. In 
vacuum, the probability that a neutrino of flavor a and 



energy E v will be detected a distance L from the source 
as a neutrino of flavor (3 is given by 



sin 2 ipjh (1) 



j,k = l 



with 



and Ajk = raj — mi, where 



ifjk ■= 1.27 A jk L/E u and A jk = to 2 - m k , 
L is measured in km, E v in GeV, and the mass eigen- 
values rxii in eV. If (anti-)neutrinos travel an appreciable 
distance through matter of sufficient density, then one 
must add to the Hamiltonian an effective potential to 
account for the (anti-)neutrino-matter interactions [13] • 
This potential yields different effective mixing angles and 
neutrino masses. For neutrinos which propagate through 
the earth over long baselines, it is crucial to take such 
matter effects into account. We do so by using the ap- 
proach developed in Ref. [HI]. We employ a two den- 
sity model of the earth: a mantle of density 4.5 gm/cm 3 
and a core of density 11.5 gm/cm 3 with radius 3486 km. 
This approach incorporate the MSW effect into a three- 
neutrino framework without the use of approximate os- 
cillation formulae. We note that it is possible for para- 
metric resonances to occur when neutrinos pass between 
regions of differing densities [19]; our treatment of mat- 
ter effects automatically insures that such resonances are 
fully incorporated. 

Atmospheric neutrino experiments are unique in that 
the baselines span several orders of magnitude making 
them sensitive to an enormous region of relevant param- 
eter space. However, the SK-atmospheric experiment is 
the most difficult to model as one must use the detected 
charged lepton to infer the direction and energy of the 
incident neutrino. A complete description of our analy- 
sis tool can be found in Ref. [2(J. Statistical errors are 
included in the chi-square function whereas systematic 
errors are accounted for by using the pull method [21 1. 



We include 43 pulls, the most important being the overall 
flux normalization. We also include [13] a simple model 
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FIG. 1: [color online] Ax 2 versus #13 in the sub-dominant ap- 
proximation. The solid [black] curve utilizes only atmospheric 
data; the dashed [blue] curve contains the atmospheric, LBL, 
and CHOOZ experiments. 

of the multi-ring events. For CHOOZ, K2K, and MINOS, 
we utilize standard analysis techniques [201 ]. 

We introduce the commonly used "sub-dominant ap- 
proximation" to provide a comparison for the full three- 
neutrino treatment. It arises from an expansion in the 
ratio of the mass-squared differences, A12/A32. The os- 
cillation probabilities are then given by 

Vee = 1 - sin 2 26<i3 sin 2 (^ 32 ) 

V efl = sin 2 6» 23 sin 2 20 13 sin 2 (ip 32 ) 

V m = 1 - 4 cos 2 013 sin 2 fe 3 (1 - cos 2 i3 sin 2 23 ) 

x sin 2 (^ 32 ) . (2) 

Additional correction terms can be added @, We 
effect this approximation by setting A21 = in our full 
three-neutrino code; in this treatment, matter effects will 
differ slightly from the approximations used by others. 

We take the bounds on the mixing angles as 0i 3 G 
[— 7r/2, 7t/2] and 812,823 G [0, 7r/2], as suggested in 
Ref. [23|. For no CP violation, this produces an allowed 
parameter space that is a single connected region. The 
(equivalent) often used bounds, 8jk G [0, 7r/2] with Dirac 
CP phase S = 0, 7r, produce two disconnected regions. 

We begin utilizing the sub-dominant approximation. 
Ref. [3] showed that atmospheric data alone restrict the 
allowed value of #13, although less so than does CHOOZ. 
We use the data from Ref. [3], which are binned more 
finely in energy than the original data Q. In Fig. [U 
we plot Ax 2 versus 6*13, varying 6*23 and A23. The solid 
[black] curve contains only atmospheric data. Our re- 
sults quantitatively reproduce those of Ref. [H which 
also used the sub-dominant approximation. Both give 
sin 2 0i3 < 0.14 (or |0 13 | < 0.38) for A X 2 < 4.6. Re- 
producing this result is a strong test for our analysis 
tool. The dashed [blue] curve in Fig. [1] represents Ax 2 
for the data set which includes the atmospheric Q , LBL 




FIG. 2: [color online] A\ 2 versus #13. The dashed [blue] curve 
includes all data and uses the sub-dominant approximation. 
The solid [black] curve incorporates the same data with a full 
three-neutrino calculation. Ax 2 for both cuves is normed by 
the minimum value of the full three-neutrino calculation. 



(K2K and MINOS [2]), and CHOOZ [| experiments. It 
is CHOOZ which restricts 6*13 much more so than the 
atmospheric data in the sub-dominant approximation. 

We now compare the sub-dominant approximation 
with the full three-neutrino calculation. In Fig. [21 we 
plot Ax 2 for the full data set (atmospheric, LBL, and 
CHOOZ) using the sub-dominant approximation, dashed 
[blue] curve, and the full three-neutrino calculation, solid 
[black] curve. We fix the solar parameters at their best 
fit values @, 0i2 = 0.58 and A 2 i = 8.0 x 10" 5 eV 2 , as 
this analysis largely decouples from these two parame- 
ters. Note that the sub-dominant results are symmetric 
about 0i3 = 0, as is manifest in Eqs. ([2]). For the full 
three-neutrino model, the asymmetry about zero is pri- 
marily due to terms in the oscillation probabilities which 
are linear in 0i3 . We have previously investigated the im- 
portance of such terms for very long baselines, L/E > 10 4 
km/GeV 0,1- Such linear terms and the interference 
between the A21 and A 32 oscillations have also been in- 
vestigated in Refs. [tj [l(|. The importance of this very 
long baseline region in the atmospheric data has also been 
noted in Ref. [24j. Terms linear in 0i3 are not negligible 
when extracting the value of this small angle. 

We next examine the contribution of the atmospheric 
data alone to Ax 2 , the [black] solid curve in Fig. [3] The 
dashed [blue] curve employs the full data set. For positive 
013, the atmospheric data are more restrictive than even 
CHOOZ. The restrictions for negative 0i3 are set by the 
CHOOZ data. Overall, we find the allowed region for 0i3 
to be asymmetric about zero, bounded from above by 
atmospheric data and bounded from below by CHOOZ. 
The final value is 8 13 = -0.07±£" at 90% confidence 
level, corresponding to Ax 2 = 6.25 for a three parameter 
analysis of this data set. In Fig. [31 we, as have others 
[25I ] . find a non-zero value for 0i3; furthermore, we find a 
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FIG. 3: [color online] Ax versus 813 using a full three- 
neutrino model. The solid [black] curve utilizes only at- 
mospheric data; the dashed [blue] curve contains the atmo- 
spheric, LBL, and CHOOZ experiments. 



statistically insignificant preference for a negative value. 

Which subset of atmospheric data results in the strict 
upper bound on #13 and the lack thereof from below? To 
answer this, we examine # 13 = ±0.15 which has A\ 2 ~ 9. 
We find that the sub-GeV fully contained events are re- 
sponsible for two-thirds of this Ax 2 . Furthermore, one- 
half of the total change in chi-squared (4.5) comes from 
the single angular bin, — 0.8 < cost? < —0.6, bin II, for 
fully contained e-like events, and the two lowest energy 
bins in which the charged lepton momentum, pi, is less 
than 400 MeV/c. This is well into the very long baseline 
region mentioned previously where we expect contribu- 
tions from terms linear in 6*13. Bin I, —1.0 < d < —0.8, 
contains neutrinos which traverse the core suppressing 
their amplitude of oscillation. 

In Fig. [4] we plot oscillation probabilities P Q , ( a(i5 i 7 1 ) 
for angular bin II and the lowest energy bin, pe < 240 
MeV/c. The solid curves use the best fit parameters, 
the dash (dot-dash) curves use 6*13 = +0.15 (—0.15). For 
"P ee , the top two [blue] curves, there is only a quadratic 
term in #i 3 which lowers V ee ■ Examining the lower three 
[red] curves, we see that V efi decreases (increases) with 
positive (negative) #13. The number of e-like events 
measured in an atmospheric experiment is related to 
R e = Vee + rVeyi, with r the ratio of the to v e flux 
at the source. From Fig. [4] effects due to terms linear in 
6*i3 combine constructively for positive #13 and destruc- 
tively for negative $13 . This is confirmed in Fig. [5] where 
we compare the effect of positive and negative #i 3 upon 
R e for neutrinos in bin II. Here, wc also discover that the 
excess of e-like events at low energies Q results in the 
strict bound on positive values of #13 in contrast to nega- 
tive values. The effect is enhanced by an MSW resonance 
near E v = 100 MeV for a mantle density of 4.5 gm/cm 3 . 
Notice that the effect is for bins where the neutrino trav- 
els through almost the entire Earth, a distance such that 
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FIG. 4: [color online] The oscillation probabilities Vee and 
Ve^ versus the inverse neutrino energy E^ 1 . The probabilities 
have been averaged over the angular bin II and folded with 
a 6% error in the energy. The solid curves use the best fit 
values of the parameters. The dash (dot-dash) curves are 
0i 3 =+0.15 (-0.15). 



the resonance is fully developed. Previously, the con- 
stanc y o f R e was found to impose an upper bound on 
I # 13 1 [26j; here we see that R e also helps constrain the 
sign of this mixing angle. 

Probing further, we may approximate the e-like events 
at Super-K for the very LBL sub-GeV data as in Ref. [§] 



R< 
1 



r sin 2 26? 2 [\ 



cot(20S)0i3-s] w?<&, (3) 



where the superscript m refers to the effective parame- 
ter values in matter and e represents the deviation from 
maximal mixing, #23 — 7r/4 + e. Note the term linear in 
613 is proportional to cot 26^. At the MSW resonance, 
8" 2 = 7r/4, and this mixing angle increases with energy 
up to 7r/2. Thus the coefficient of the #i 3 term in Eq. ([3|) 
is negative above the resonant energy so that a negative 
6*i3 produces an excess of e-like events. The bounding 
of 6*i3 from above by the atmospheric data depends criti- 
cally on incorporating the full MSW effect. For these low 
energy neutrinos, one has r ~ 2; the dependence of our 
results upon this ratio is not severe. More significant is 
the dependence upon e; the octant of 6*23 can conspire to 
enhance or suppress the effect of #13 upon an excess of 
e-like events. 

In this new era of precision neutrino experiments, small 
effects, such as those arising from #13, require a careful 



treatment. Future reactor experiments [111 . |12|, [13J are 



sensitive to # 2 3 and thus can determine the magnitude of 
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#13 , but not its sign. The long baseline experiments 
1 151 ] will contain small effects that are linear in #i 3 , while 
an upgraded Super-K will produce additional data in the 
region most sensitive to effects linear in #13. How these 
different data interplay with each other in determining 



4 




1/E (GeV 1 ) 



FIG. 5: [color online] The measured quantity R e versus the in- 
verse neutrino energy E^ 1 for angular bin II. The solid [black] 
curve utilizes the best fit parameters, the dashed [blue] curve 
813 = +0.15, the dot-dashed [red] curve #13 = —0.15 



this mixing angle will be most interesting. 

We find that present atmospheric data restrict the 
value of #13 from above, while the limit from below 
remains as determined by CHOOZ, and that #13 = 
—0.07^0 11, assuming no CP violation. It is important 
to realize first that #13 can be negative [23[ and, second, 
that linear effects lead to asymmetric errors. Our analy- 
sis requires the use of the more finely binned atmospheric 
data, Ref. 0] , the use of the full three- neutrino oscillation 
probabilities, and inclusion of the full MSW effect. 
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